Catching the zitterbewegung 
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We demonstrate how the zitterbewegung charge oscillations can be detected through a charge con- 
ductance measurement in a three-terminal junction. By tuning the spin-orbit interaction strength or 
an external magnetic field the zitterbewegung period can be modulated, translating into complemen- 
tary conductance oscillations in the two outgoing leads of the junction. The proposed experimental 
setup is within the reach of demonstrated technology and material parameters, and enables the 
observation of the so far elusive zitterbewegung phenomenon. 
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I. INTRODUCTION 

The quantum relativistic theory of the electron, ex- 
plaining the "peculiar two-valuedness" or spin angu- 
lar momentum, was put forward in 1928 by Dirac. 
While studying free particle solutions to Dirac's theory, 
Schrodinger found that the time dependence of the po- 
sition operator of an electron traveling through vacuum 
displays an oscillatory motion around some mean posi- 
tion, named the zitterbewegung, arising due to an inter- 
ference between the positive and negative energy states 
of the Dirac electron. The angular frequency of the os- 
cillation is oj = 2mc 2 /h and the amplitude is twice the 
Compton radius A = c/lu = 2r c = h/mc, implying an 
angular momentum s = mcr c — h/2. This equality 
of the zitterbewegung amplitude and the intrinsic spin 
angular momentum of the electron has led to the sug- 
gestion that indeed the zitterbewegung is the fundamen- 
tal mechanism behind the intrinsic angular momentum. 
Hestenesi extended this idea further by suggesting that 
the zitterbewegung is a truly fundamental property of the 
entire quantum reality responsible for the appearance of 
the complex phase. Establishing whether the zitterbewe- 
gung is indeed a real physical observable is therefore of 
great importance for a deep understanding of the quan- 
tum reality. However, due to the extraordinary frequen- 
cies and amplitudes set by the electron-positron energy 
gap, a direct observation of the electron zitterbewegung 
in vacuum is beyond the reach of standard measurement 
techniques. 

In the non-relativistic limit the intrinsic angular mo- 
mentum of the electron is a two-component spinor prop- 
erty. Here, the positive-negative energy interactions of 
the four-component Dirac equation are folded into spin- 
dependent terms in the Hamiltonian. For an electron 
with non-zero momentum, p, in an electric field, E, this 
gives rise to a coupling of the spin and spatial degrees of 
freedom through the spin-orbit interaction (SOI) Hamil- 
tonian ifsoi = (efi/4m 2 c 2 )(T • E x p.— The SOI gives rise 
to interference properties between the two spin states in 
a way analogous to the interference between the positive- 



and negative-energy wave components of the relativistic 
Dirac particle in vacuum. The prefactor contains the in- 
verse of the electron-positron creation energy (« 1 MeV) 
leading to the difficulty in the observation of effects due 
to the SOI for a free particle in vacuum. In a semicon- 
ducting crystal, however, the relevant energy scale is the 
conduction- valence band gap which is about six orders of 
magnitude smaller than the electron-positron gap, mak- 
ing SOI effects of strong influence on the spin dynam- 
ics. One thus heuristically expects the zitterbewegung 
to be of much larger amplitude and lower frequency in 
crystalline materials^^ Here, the electric field which 
induces the SOI is naturally supplied by the atomic po- 
tentials as well as interfacial potentials between different 
crystalline materials. The SOI derived from the atomic 
potentials is called the Dresselhaus SOI and is naturally 
a fixed material property. The SOI derived from the in- 
terfacial field is called the Rashba SOI and can be tuned 
by, e.g., a voltage applied to a surface gate^iS tilting 
the band structure of a quantum well. It was indepen- 
dently shown in Refs. [H, d, and [h] that zitterbewegung 
charge oscillations on the scale of 100 nm arise in a semi- 
conducting waveguide with a spin-polarized electron in- 
jection in the presence of Rashba SOI. Large spin po- 
larization has been theoretically predicte d 11 : 12 and ex- 
perimentally achieved in three-terminal semiconducting 
structures M The spin separation in these structures rely 
on the spin-Hall effect induced by the SOI. It was shown 
in Ref. |9j that the zitterbewegung and the spin Hall effect 
are intrinsically connected and only differs in injection 
and detection conditions. It should therefore be possible 
to measure the zitterbewegung in these structures. 

Calculations of the average position of a free (un- 
bound) Dirac particle^ as well as wave packets in 
graphene layers and carbon nanotubes^ showed that 
the zitterbewegung oscillations have a transient character 
with a decay time of ~ 10~ 12 — 10~ 15 s. This transient 
decay occurs due to a spreading of the wave packet, as 
well as separation of the subpackets corresponding to the 
positive and negative energy solutions. However, Lochia 
pointed out that in a confined system this spreading is 
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FIG. 1: Schematic description of the considered three- 
terminal junction structure. Arrows indicate schematically 
the spin-polarization evolution and the quasi-classical zitter- 
bewegung path of electrons. 



restricted and the transient character would therefore 
not appear for bound states. An attempt to observe 
the zitterbewegung should therefore be directed to quan- 
tum confined systems where periodic oscillations can be 
sustained. The approach of studying wave-packet solu- 
tions might not be very well suited for the steady-state 
electron flux in the crystal. Indeed, the electron flux 
is not a localized "particle" but rather a quasi-particlc 
excitation represented by an extended state. A proper 
analysis of these extended state solutions are thus called 
for. Previous suggestions for measurements in quantum 
wires include direct scanning gat e 16 ' 17 ' 18 ' 19 or Kerr ro- 
tation measurement s 20 ' 21 of the zitterbewegung oscilla- 
tions. This would require resolutions of the order of tens 
of nanometers, which is on the cutting edge of scanning 
gate techniques and far beyond the capabilities of the 
Kerr rotation method. 

In this paper, we will analytically derive through sec- 
ond order perturbation theory an expression for the zit- 
terbewegung charge oscillation patterns for steady-state 
quantum transport in a quantum wire. We also perform 
exact numerical calculations of the transport properties 
for a three-terminal junction. We will show how a simple 
charge conductance measurement can be performed to 
measure the zitterbewegung in such a junction defined in 
a semiconducting material. The proposed experimental 
setup, illustrated in Fig. [TJ employs standard nano-scale 
lithographic techniques as well as experimentally demon- 
strated material parameters. 



the Hamiltonian (with H = 1) 



Ho + H s , 



(1) 



where Hs contains all the spin-dependent terms. We will 
choose the spin quantization axis to be along the mag- 
netic field direction. The scattering state at a given en- 
ergy in the quantum wire can be written as an expansion 
of the eigenstates of the wire 



(2) 



The analytical solutions for the eigenstates of the spin- 
orbit coupled quantum wire are not know. We therefore 
construct the states from second-order perturbation the- 
ory where the unperturbed states are the eigenstates of 
H : 
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The unperturbed eigenstates are spin degenerate. How- 
ever, the total Hamiltonian is symmetric under the 
cr y Ry operation (where R y is the reflection operator 
in the transverse direction). This operation defines 
a pseudo-spin (generalized parity) which block diago- 
nalizes the Hamiltonian. We can therefore perform 
the calculations separately in the pseudo-spin spaces 
s+ = {|1T^),|2|(°>),|3T (0) ),|4|(°)),...} and a~ = 

{|U(°)),|2T^),|3l(°)),|4tW),...}. 

The general form of the first-order correction to the 
pseudo-spin states read 
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where we have introduced the shorthand notation 
Cicrja' = (i,a^\aa x p y + a y (g*^ B B y /2 - ap x )\j,a'^) € 
R. The general form of the second-order correction to 
the s 1 * 1 pseudo-spin states read 



II. 



ZITTERBEWEGUNG OF THE ELECTRON £ s ±(2) ) = 

FLUX IN A QUANTUM WIRE 



We start by considering a uniform quantum wire (the 
left part of Fig. [I]). The considered quantum wire is de- 
fined by a hard-wall confinement potential and is sub- 
jected to the Rashba SOI and a magnetic field B y ap- 
plied along the y direction. The system is described by 
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Writing out explicitly the first-order correction to the 
lowest s + pseudo-spin state we have 
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where we have truncated the converging series and set 
Ei = ir 2 /2m*w 2 . The second term in Eq. 0, c 3 jif , van- 
ishes due to the orthogonality of the unperturbed states. 
Similarly, for the lowest s~ pseudo-spin state, the first 
order correction reads 
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where in the last equality we have used the fact that 
Ciaja' — ~ c iaja'i with the bar indicating a spin flip. 
We can then put both first order corrections to the two 
pseudo-spin states on a short form as 
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The second order correction to the lowest s + pseudo-spin 
state reads 
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A similar expression can also be derived for the second or- 
der correction to the lowest s~ pseudo-spin state. Here, 
due to the second order nature, the coefficients are in- 
variant under the spin flip and we can thus put the two 
second-order corrections on the short form as 

„(2), 



|l,s+( 2 )) 



4 2) |2|(°>) 



■4 2) |3T (0) ) 



(0)\ 



|l s -( 2 ') « C ( 2) |2T (0) )+4 2) |3| (0) )+4 2) |4T (0) >. 

(11) 

We want to study the zitterbewegung pattern for trans- 
port in the two lowest eigenchannels of the wire and, 
therefore, calculate the expectation value of the trans- 
verse coordinate, 
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But since the spatial coordinate operator y is of odd par- 
ity we have (1, s + \y\l, s + ) = (l,s~\y\l,s~) = 0. Thus, 
Eq. (|12p implies that there is no uniform displacement 
of the wave function in the y-direction derived from the 
spin-orbit interaction or the in-plane magnetic field B y 
(this is true to all orders). We are hence left with calcu- 
lating the term 
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The first term is obviously zero. From Eq. §§§ we can 
easily see that the second and third terms are mutually 
cancelling. Furthermore, from Eq. © we see that the 
first order corrections to the lowest s + states only con- 
tains unperturbed spin-down components and vice versa 
for the corrections to the lowest s~ states. This means 
that the fourth term is identically zero. From Eq. pip 
we see that the fifth and sixth terms in Eq. (|13[) combine 
to a second-order non-zero contribution. The seventh, 
eighth, and ninth terms are third- and fourth-order con- 
tributions and will hence be discarded. We are then left 
with the result 
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The second term can safely be neglected since it repre- 
sents a correction of the order c\ /c y 2 ' = 4/5 . 

The coefficients a a in Eq. (|12[) are determined through 
the injection boundary condition. We consider the con- 
dition that electrons are injected into the wire from a 
non-interacting semi-infinite lead. We further assume 
that the spin-orbit interaction is adiabatically increased 
towards the finite value in the wire. If the injected elec- 
trons are polarized in the out-of plane z direction we 
have a z ls+ = a^ s _ = 1/V2- Similarly, for electrons in- 
jected with a polarization in the longitudinal x direction 
we have a? + — 1/V2 and a? _ = i/^/2 and for electrons 
polarized in the transverse y direction we have a^ s+ = 1 
and aO. _ = 0. This gives the zitterbewegung pattern 
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Our derivation of Eq. (fT5|) has shown that the steady- 
state zitterbewegung pattern for the electron flux in a 
quantum wire is a second order process in the spin-orbit 
interaction. It arises through interference between the 
spin-split eigenstates at a finite k vector. The amplitude 
of the oscillations depends on the SOI strength, the mag- 
netic field applied, and the quantization energy of the 
wire. Furthermore, it is periodic with a period propor- 
tional to the inverse of the k vector difference between 
the split states, 1/Ak. In order to find this period we 
calculate the eigenenergy corrections to second order 

rp „ 7T 2 fc2 , 9*^bB v a 2 ™* 

which gives us 

ki s + ss am* 
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Similarly, for the lowest s~ pseudo-spin state, we have 
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V w 

The linear response transport is taking place at the Fermi 
energy Ep. Setting both pseudo-spin energies to Ep and 
considering the positive branch we have the wave-vector 
difference 

fc ls + — k\ s - = Afc « 2am* 
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V ur 
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In the limit of 2m*E F - ^ > |2(am*) 2 ± m*g*/i B B y \ 
we can expand the square roots to first order giving 

U*m*(2a- - €m£i ) (20) 
\ yj2m*E F - n 2 /w 2 J 

The zitterbewegung period can hence be tuned by chang- 
ing the SOI strength as well as the external in-plane mag- 
netic field. The magnetic field alone causes a rotation of 
the spin angular momentum (for x or z polarized injec- 
tions) but does not generate a spatial charge oscillation. 
When the spin and orbital degrees are coupled through 
the SOI, the rotation of the spin angular momentum cou- 
ples into a spatial charge oscillation. Quantum mechani- 
cally this manifests itself in interference between the spin- 
split states at the Fermi energy. The degree of splitting 
is tuned by the SOI strength and the magnetic field, cre- 
ating a modulation in the oscillation period. 



III. ZITTERBEWEGUNG IN A 
THREE-TERMINAL JUNCTION: EXACT 
NUMERICAL SOLUTIONS 

In order to devise an experiment to measure this ef- 
fect we consider the three-terminal junction structure, 
illustrated in Fig. [U to be defined in a two-dimensional 
electron gas formed inside a semiconductor heterostruc- 
ture. Manufacturing of such quantum structures can 
routinely be made by, e.g., state of the art lithographic 
techniques i 13 ' 22 i 23 ' 24 On top of the structure, a metal 
gate is placed, enabling tunability of the electric field 
in the heterostructure and, hence, the Rashba SOI 
strength. Tunability of the Rashba SOI strength in the 
ranges of 6-9 meV nm and of 1.5-5 meV nm has been 
demonstrated in InGaAs quantum wells. 6,7 The junc- 
tion is injected with spin-polarized electrons from a spin- 
filter devicei 11 ' 12 i 13 ' 25 i 26 ' 27 i 28 ' 29 The waveguide presents 
a quantum confined structure in the transverse direc- 
tion eliminating the detrimental effects of zitterbewe- 
gung transients expected in extended systemSi 14 i 15 The 
results of Refs. d, and [1(3 as well as Eq. (fT5|) show that 
the injected carriers will exhibit a transverse oscillatory 
charge pattern in the waveguide. Intuitively, depending 
on the period of oscillations the carriers will then be cap- 
tured with unequal probabilities in the upper and lower 
branches of the junction, as indicated in Fig. [T] Now, 
according to Eq. (|20[) the period of oscillations can be 
tuned through a or B y and we therefore expect a peri- 
odic modulation of scattering into the upper and lower 
branches, respectively, as we tune the SOI strength or an 
external magnetic field. By measuring the conductances 
between the central and the upper branches and between 
the central and lower branches, respectively, we should be 
able to detect the zitterbewegung oscillations as a com- 
plimentary oscillating conductance signal. To study this 
we calculate the spin-resolved quantum transport prop- 
erties in such a three-terminal structure by adopting an 
exact numerical scattering-matrix method developed in 
Refs. HHIHIM 

The SOI is non-zero and spatially uniform in the cen- 
tral region of the structure shown in Fig. [T] The perfect 
leads (not shown in Fig. []} attached to the branches are 
assumed to have vanishing SOI. Since we shall consider 
charge transport, this assumption will only induce slight 
scattering at the lead-branch interface and will hence not 
change the physics we are studying here. Figure [2] shows 
the numerically calculated charge conductance into the 
two outgoing leads (solid and dashed lines, respectively) 
of the three-terminal structure at B y = as a func- 
tion of the SOI strength. The different solid lines rep- 
resents varying degree of polarization of the injected car- 
riers ranging from a full spin polarization p = 100% to 
a weak spin polarization p — 20%. As the SOI strength 
is modulated the charge conductance displays an oscil- 
lating behavior. The oscillations of the charge conduc- 
tances into the upper and lower leads are exactly out of 
phase and increase in amplitude with increasing the SOI 
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FIG. 2: Charge conductance at B y = into the upper (dashed 
line) and lower (solid line) leads as a function of the SOI 
strength, for carriers injected into the central lead with differ- 
ent degrees of spin polarization along the z direction. The zit- 
terbewegung oscillations are detected as two out-of-phase con- 
ductance oscillations. Spin-orbit interaction strength is given 
in terms of the width w and the lowest transverse eigenen- 
ergy E\ = tt 2 /2m*w 2 of the leads. The Fermi energy is set 
at E F = 3 Ei. 




FIG. 3: Charge distributions for an electron injection with 
an initial spin-polarization along the z-direction in the three- 
terminal junction with the SOI strength of (a) a — 0.129 wE\ 
and (b) a — 0.199 wE\. The change in the SOI strength 
results in a modulation in the zitterbewegung period, leading 
to a capture of electrons in the upper or lower branch. 



strength. Going from a ~ 0.13 wE\ to a s» 0.2 wE\ 
half a period of oscillation can be achieved. For an In- 
GaAs waveguide of width w — 150 nm, this corresponds 
to 6.5 meV nm < a < 10 meV nm, translating to a gate 
voltage sweep of AV g m 7 V (Ref. 0), within the reach 
of the tunability demonstrated in Ref. H. 

To demonstrate that these alternating conductance os- 
cillations arise due to the modulation of the zitterbewe- 
gung period, we plot in Fig. [3] the charge probability dis- 
tribution of the electron wave in the three-terminal junc- 
tion structure. In Fig.^a) the calculation is made for the 
SOI strength a = 0.129 wE±, i.e., at an extreme point 
of the conductance to the upper lead. As the initially 




-10 -5 5 10 

k [1/w] 

FIG. 4: (a) Charge conductance into the upper (dashed line) 
and lower (solid line) leads as a function of the magnetic field, 
for carriers with initial spin-polarization along the z-direction 
injected into the central lead at the Fermi energy Ef = J>E\. 
The zitterbewegung oscillations are detected as two out-of- 
phase conductance oscillations. At B y ~ 0.88 Ei/g*fj,B 
(marked by arrow) the SOI and Zeeman fields are mutually 
cancelling, creating equal conductances, (b) Band structure 
of a uniform wire for B = 0.88 g* fiB /Ei. At the Fermi energy 
Ef = 3 Ei (dashed line), the k vector difference of the two 
lowest forward-propagating states is Ak = due to the can- 
celling SOI and Zeeman fields. In both panels a = 0.1 wEi. 



spin-polarized electrons travel along the central branch 
the SOI induces zitterbewegung oscillations in the charge 
distribution. When the electrons reach the junction, the 
central divider captures the oscillations into the upper 
branch, resulting in a high conductance to that lead. 
Conversely, in Fig. [3]Jb) the calculation is made for the 
SOI strength a = 0.199 wE\, i.e., at the half-period shift 
of the extreme point of the conductance to the upper 
lead. Here the zitterbewegung charge oscillation period 
is increased and the central divider captures the oscilla- 
tions into the lower branch. 

The zitterbewegung period can also be tuned by the 
in-plane magnetic field. In Fig. 0|a) we show the cal- 
culated conductances into the two outgoing branches 
as a function of the magnetic field at the fixed SOI 
strength a = 0.1 wE\. At positive fields the conduc- 
tance show regular oscillations with changing B y . At 
By rj 0.88 E\/g*HB (marked by an arrow) we can see 
that the conductance into both outgoing branches are 
equal (two extreme points). We can understand this 
from Eq. (|15[) which says that the zitterbewegung am- 
plitude is proportional to the vector addition of the SOI 
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FIG. 5: Same as Fig. [2] except that the Fermi energy is set 
at Ef = 12 Ei, corresponding to a multi-channel electron 
transport through the three-terminal junction structure. 

and Zeeman fields. At the point 2ak — g* [IbB v the 
two fields are mutually cancelling and the zitterbewe- 
gung patterns disappear. Taking the k vector to be the 
unperturbed value at the Fermi energy, this occurs at 
By = 2air^/2/wg* hb ~ 0.89 E 1 /g*fi B - This explains the 
equal conductance into both outgoing branches. Since 
the fields are cancelling we also expect the zitterbewe- 
gung period, given by 2tt/ Ak, to be infinite at this point. 
By inspecting the band structure for a quantum wire 
[found by exact numerical diagonalization of Eq. I|T|)] 
shown in Fig. [31(b) we can see that at the Fermi en- 
ergy (marked by a dashed line) the two forward prop- 
agating states with positive k number have Ak = at 
By 0.88 E\/g*fiB, which indeed implies an infinite 
zitterbewegung period. The suppression of the zitter- 
bewegung oscillations at this point is hence related to 
the mutual cancellation of the fields and the subsequent 
restoration of the time-reversal symmetry and recovery 
of spin as a good quantum number in the system^ In 
the negative magnetic field domain the SOI and Zeeman 
fields are mutually constructive and there is hence no 
such degeneracy point. 

The above calculations were in the single-channel 



transport regime. For a multi-channel transport 
the zitterbewegung patterns becomes increasingly 
complicated. 9 However, from the calculated conduc- 
tances in the multi-channel transport regime (Fig. [5]) , we 
see that the tell-tale complementary oscillations akin to 
the zitterbewegung can nevertheless still be clearly dis- 
cerned. 



IV. CONCLUSION 

We analytically derived the zitterbewegung charge pat- 
terns for the steady-state electron flux in a quantum wire. 
The zitterbewegung is a second order process in the spin- 
orbit interaction and its spatial period and amplitude can 
be shifted by tuning the spin-orbit interaction strength or 
an external magnetic field. This was confirmed by exact 
numerical calculations of the charge distribution and the 
conductance. By performing a charge conductance mea- 
surement in a three-terminal semiconductor nanostruc- 
ture it is possible to detect the zitterbewegung charge 
oscillations, that has eluded researchers for almost 80 
years, as a complimentary conductance oscillation in the 
outgoing arms. This oscillation should be detectable even 
for injected carriers with spin-polarization as low as 20%. 
Due to the relative simplicity of the set up and the read- 
ily available enabling technologies we expect that such a 
measurement could be performed immediately. Our pro- 
posal has a great advantage to, e.g., a two-dimensional 
electron system defined in a quantum well, which suffers 
from time decay of the zitterbewegung, and can sustain 
periodic oscillations due to the quantum confinement. 
By performing the suggested experiment the question re- 
garding the reality of the zitterbewegung could be settled, 
and thus bring new light on the very underlying nature 
of quantum reality. 
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